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FOREWORD 

This is the final reporting on a one man-year study of techniques 
to analyze aerodynamic noise sources in a free field and in the reverber- 
ant environment of the 40 x 80 foot wind tunnel at Ames Research Center 
(ARC). The project was funded by ARC-NASA, under contract No. NAS 2-6872. 

The goal of this project, in a broad sense, was to identify and 
define techniques for studying aerodynamic noise generating mechanisms 
without perturbing the flow. Identification and definition of noise 
sources entails localization of such sources and determination of their 
characteristics (spectrum, strength, coherence, phase, etc.) under various 
flow regimes. Both analytical and experimental investigations were to be 
conducted. The experimental task was to be done in the 40 x 80 foot wind 
tunnel at ARC; consequently, the adopted technique of measurement was to 
be capable of discriminating between the true aerodynamic noise sources 
and the reverberant or background noise in that tunnel. 


vii 
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SUMMARY 

This report covers: (a) phase I, the theoretical analysis of the 

various diagnostic techniques to define and identify noise sources in 
free field and, more importantly in ARC 40 x 80 foot wind tunnel which is 
quite reverberant, and (b) phase II, the design and testing of an acoustic 
antenna with an electronic steering control. Such a scheme for noise 
diagnosis was agreed upon because of its relatively low cost, short time 
needed for its design and construction and its promised good performance. 

Two models (909 and 909B) of the acoustic antenna electronic steering 
control were built. The second model has higher directivity (10° half 
angular resolution at 6 dB) as well as higher frequency selectivity. The 
antenna utilizes two crossed linear arrays of microphones, (5 microphones 
in each array) ; each array can be steered independently from 0° to 180° . 

The other important feature of the design is that the steering is indepen-^ 
dent of the acoustic frequency. The output in the first model is the sum 
of the outputs of the two arrays; the second model measures the product 
of the two arrays. 

In the course of study, the theory of the acoustic characteristics 
and reverberation modes of wind tunnels was developed to show that a wind 
tunnel is quasi-reverberant , having a characteristic length l of the order 
of magnitude of a wavelength. For a point source, on the axis of the tunnel, 
approximately plane waves propagate beyond and higher order modes attenu- 
ate rapidly within the volume specified by & (Sec. 1.2.1). Extended and/or 
non-axial sources result in high order modes propagating down the tunnel. 
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In all cases, however, there is a quasi-reverberant field within the 
volume characterized by the length &. Thus, estimating the true strength 
of a sound source radiating in the reverberant environment of wind tunnels 
requires measurement in extreme proximity of the source and hence highly 
directional schemes. 

In addition, calculation was initiated to obtain the reflection 
caustics of various orders in the tunnel; however, the work was abandoned 
in favour of pursuing the wave theory to obtain the mode of reverberation 
and propagation. 

Techniques (mainly optical) to correlate the results of the 
acoustical antenna noise measurements with the density fluctuations in 
a jet flow are pointed out for future implementation. 

The theoretical phase of the study is detailed in parts I, II, III 
and Appendix A of the report. The discussion on the acoustical antenna 
and its electronic steering control is covered in part IV and Appendix B. 
Conclusion and recommendations for further extension of the project are 
expounded in part V. Two bibliographies appear following part V and 
Appendix A. 

Units used throughout the report are the International System of 
Units and the customary units. 


Valuable technical support in one aspect or another of this study was 
contributed by K. Bracewell, D. Bershader and P. Soderman, to whom the 
authors are indeed grateful. 



1 . 


EL-SUM CONSULTANTS 


NAS 2-6872 


GENERAL INTRODUCTION 

In spite of recent advances in noise control engineering, it is 
good engineering practice^^ to perform full size noise tests on aircraft 
engines and/or components. These tests will allow the engineer to 
introduce final improving adjustments to the design, since our scaling laws 
for noise generation are not that precise. 

Now, it is important to obtain the noise characteristics of these 
engines (or components) not only when the engine is idling on the ground, 
but also under flight condition. Whereas, engine testing in a wind tunnel 
provides satisfactory information about the aerodynamic behavior of the 
engine, unfortunately for the acoustical engineer conventional noise 
measurements in a wind tunnel yield data which are difficult to unfold 
for the purpose of extracting the free field noise characteristics of 
the engine being tested. When the engine is in a wind tunnel, the engine 
noise radiation pattern as well as its spectrum are both distorted. At 
low Mach numbers the acoustic power emitted by the engine, however, is 
not affected by its being placed inside the wind tunnel, if the tunnel 
test section is large enough so as not to affect the aerodynamic per- 
formance of the engine. 

One of the main reasons which makes it difficult to obtain the 
free field noise radiation pattern from wind tunnel measurements is the 
interfering presence of the walls. Indeed, in addition to the noise 
signal directly radiated from the engine under test, a detector would 
pick up signals reflected from the walls. Since the walls of the wind 
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tunnel are fairly good reflectors, this reflector or reverberant signal 
can be comparable to the directly radiated signal. It is important to 
point out that had one used an acoustic wattmeter of highly directional 
characteristics, the directivity pattern of the source (i.e., engine 
under test) could have been obtained. Unfortunately, at the present, no 
acoustic wattmeters with such desirable highly directional properties 
are available, so that one has to resort to other means such as are des- 
cribed later in this report to measure the direct sound radiated. Another 
reason responsible for complicating the interpretation of wind tunnel 

noise testing is the noise created by the fans driving the air inside 

( 2 ) 

the tunnel. This broadband noise not only represents a background 
noise, but it can interact with the source thus enhancing the noise 
signal sensed by the noise detectors (see pp. A-ll , 12). 

In order to achieve the main goals of this investigation, viz., 
the development of opto/acoustical diagnostic tools to study noise 
generation without perturbing the flow, first we have to understand the 
interpretation of the measurements made in the tunnel. Accordingly, we 
propose to select an engine whose aero-acoustic characteristics are 
fairly well understood. The general noise features of such an engine 
in the tunnel are then outlined. Now a turbojet run cold (say by moving 
the compressor blades by an external electric motor) meets many of our 
requirements : 

(1) The acoustic characteristics of the resulting subsonic jet are 

(3) 

well understood. 

(2) The noise produced by the blades of the rotating machinery 
of the turbojet is rich in harmonics and appears as discrete frequencies 
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superimposed over the wideband spectrum of the jet* The mechanism for 

this noise is that of the general class of axial flow compressor noise 

(4) 

whose theory has been worked out in details* 

(3) A turbojet rather than a conventional jet flow is proposed 
because the discrete blade passage frequency of the turbojet can be used 
as a calibration signal for the diagnostic techniques we wish to use and 
will be explained later in this report* 

A fundamental question we wish to emphasize here is that a highly 
directional microphone which can be used to discriminate between the 
direct sound radiated and the reverberant noise will allow us to evaluate 
the total acoustic power generated by the turbojet. As mentioned pre- 
viously this acoustic power is unaffected by the presence of the walls 
of the tunnel* 

Another important point is that the turbulent fluctuations in the 
jet also are unaffected by the presence of the walls. Consequently, the 
measurement of the physical properties of these fluctuations (mean square 
intensity, correlation, etc.) through an optical diagnostic method will 
lead us to an estimate of the equivalent strength and spatial extension 
of the acoustic sources whose radiated power is (presumably) known. 

In part I of this report the general aero-acoustic features of 
the turbojet are described. In particular, we will summarize the scaling 
laws for noise production and the effect of the walls of the tunnel. In 
the second part of the report, we point out a suggested procedure for 
differentiating between the direct and reverberant and background noise. 
In the third part we discuss optical methods that will allow us to map 
the fluctuating (turbulent or coherent) part of the flow. Part IV deals 
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with the design characteristics of high directivity acoustic antenna 
and its electronic steering control. Two models of this antenna were 
built, tested and evaluated. The results of the evaluation are included 
in this section. Part V enumerates our recommendation of areas that 
needed further follow up. 

To make this report thorough and complete, we reproduced earlier 
suggestions and technical analysis of this project in Appendix A. 
Appendix B is a copy of the operating manual of the acoustical antenna 
and its electronic steering control. 
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PART I 

ACOUSTIC CHARACTERISTICS OF TURBOJET 
AS A NOISE SOURCE 

1.1 General Remarks for Free Field Sound Generation 

The noise generated by a turbojet engine is contributed by a 
number of processes, the two salient ones being the turbulent noise in the 
jet flow and that due to the blade-passage tones (from turbine, from 
compressor, etc.). The presence of any afterburner adds to the noise 
generation, through additional combustion noise. This latter contribution 
will be neglected here. 

1.1.1 Turbulent Jet Noise 

(5) 

Extensive measurements on turbojet noise have indicated that 
the jet noise has the following characteristics:^ 

(i) the efficiency of conversion of mechanical to acoustic energy 
is a slowly varying function of the Reynolds number at low mach numbers 
and follows an inverse fifth power law of the mach number. Actually, one 
has 

where P is the acoustic power > V the exit velocity * p Q> C q the density 
and velocity of sound at the ambient temperature and d the diameter of 

-A 

the jet nozzle* The constant K is typically of the order of 10 

An empirical relation 1 "^ which has been obtained from over two 
decades of data collection on turbojet noise yields the total acoustic 
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power as a function of a parameter £ given by 


where 



T is the exit temperature in degrees Rankine 
V is the average exit velocity in ft /sec 
p is the exhaust pressure in atmospheres 
d is the diameter of the nozzle in inches. 


( 2 ) 


(3) 


The dependence of the acoustic power P on £ is shown in graphical form 
in Fig. 1. 

(ii) The sound spectrum is very broad and has an ill-defined peak. 
The peak frequency varies as V/d. The spectral density falls off on 

( 8 ) 

both sides of the peak and a generalized spectrum is shown in Fig. 2. ' 

(iii) The polar distribution for the intensity of the sound 

radiated follows a "butterfly” type configuration in which almost all the 

sound radiated is encompassed in an acute angle with the jet. The angle 

for the directional maximum is apparently a function of the frequency 

and of the mach number. Defining a Strouhal number fd/V, Lassiter and 

( 9 } 

Hubbard state that this angle varies from 30° to 40° for fd/V values 
of 0.4 to 1.5. For lower values of this ratio, the angle slowly decreases 
to 15° for a corresponding ratio of 0.2. 

(iv) The mechanism for the generation of sound by turbulence 
becomes considerably more efficient in the presence of a mean shear. 

This increase in efficiency can be explained by the fact that the 
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Acoustic power level as a function of parameter £•> 
(Ref. 7). 
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Fig. 2. Normalized sound spectrum for a cold air jet (Ref. 8). 


momentum flux is considerably increased when it takes place in the 


direction of a large velocity gradient. Lighthill, using the result of 
Proudman^^ was able to show that the order of magnitude of the fluctua- 
tions of the gradient which contribute most to the rate of shear strain 

(~\ in / 

and hence to the sound production is about 81 v J / L, where L is the 

scale of the turbulence. Proudman had demonstrated that for isotropic 

turbulence the large nondissipating eddies contribute most to the sound 

production while the energy dissipating eddies of size much less than L 

do not radiate at all. In the presence of a mean shear the intermediate 

eddies become the dominant contributors to the sound production. 

( 12 ) 

Ribner was able to find, by combining Lighthill and Proudman 


theories to regions of similar flows for a subsonic jet, the spatial 
distribution of the equivalent noise sources along the jet. His calcula- 
tions referred to the noise power emitted by a slice of jet indicate that 
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this power is essentially constant with distance (x) in the initial 
mixing region of the jet, then further downstream about 8 to 10 diameters 
from the nozzle the power falls off very rapidly with an x ^ power law. 

In this region of the jet the turbulence is fully developed (Fig. 3). 

The bulk of the noise is hence produced in this near region to the nozzle. 
These very important calculations are in substantial agreement with 
experimental observations. 

1.1.2 Blade Passage Noise 

It is clear that any detailed prediction of the noise generated 

by the turbine and compressor of a turbojet would require an accurate 

knowledge of the configuration of these components (i.e., profile of 

blades, aerodynamic performance, etc.). What we will do instead is to 

summarize the important features of the noise production in terms of 

design parameters. This point of view incorporates the results of 

intensive research on this subject. One of the investigations we have 
(13) 

used extensively in fact includes over 100 references on blade 

passage noise. The characteristic features of this noise are: 

(i) The predominant radiation mechanism is through the fluctuating 
forces on the blades. The stator and rotor blades are equally important 
to the noise production. The coherent noise radiated is not only at the 
blades passage frequency but is rich in harmonics. Increasing the rotor- 
stator separation will preferentially reduce the higher harmonic radiation 
by the stator, while increase of stator-rotor separation will reduce all 
harmonics of the noise generated by the rotor. 

(ii) The overall power and directionality of the noise radiated 

are determined by the source parameters, the rotational speed and blade/vane 
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numbers. From a careful consideration of the fluctuating lift on the 

/ ION 

blades, Lawson v * finds that the acoustic power radiated at the funda- 
mental frequency of the blade passage frequency is given approximately by 


<>■, s BA i; B i 

128 A a 2 
e o 


(4) 


where 


is the total area of the blades 

A is the effective area based on an effective radius R, 
e 

2 

i, e. , A = ttR 
e 

is the aerodynamic drag coefficient for the blade 

is the relative velocity at some effective radius 

p ,a are the density and sound velocity respectively for 
o o 

ambient conditions and 
B is the number of blades. 


Since it is more conventional to refer the noise production as a function 
of the mechanical tip speed for the blades through the relation 


U 


1 



(5) 


where R^, is the tip radius, and when a 45° angle between and is 
assumed, one finds that Eq. (4) reduces to 


W - 


P C n B 
D 


128 a“ 




( 6 ) 


For the above expression, = D/2. It thus becomes apparent that the 
noise generation depends on the fifth power of the velocity and on the 


area of the blades. The above formula is in good agreement with 
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experimental observation. Indeed when it is compared to empirically 

(14) 

derived formulas such as that of Beranek which is 

PWL (re lCf 13 watts) = 135 + 20 log 1Q HP - 10 log 1Q q (7) 

one finds that Eqs. (6) and (7) are equivalent since 

2 3 

HP = horsepower — D V 

3 2 

and q = compressor discharge in ft /sec — D V, 

Using reasonable values for the parameters such as 

V A e ' °- 5 

K/H,, - 0.8 

C D ' °' 2 

one finds that Eqs. (6) and (7) have the same dependence on D and V. 
Equation (6) underestimates the power by a few dBs, This is expected 
since Eq. (6) is the power radiated for the fundamental only . When the 
contribution of the harmonics are added the agreement becomes very close. 

(iii) An important property (to our study) of the blade passage 
tone is the directivity of the radiation pattern. This pattern is readily 
obtained from the fluctuating rotating forces on the blades. The equiva- 
lent sources turn out to be equivalent to a ring source. Consequently 
appreciable side lobes will appear whenever the ring size becomes com- 
parable to the (Doppler-shifted) wave length. Because of the presence 
of the air flow, more sound is radiated forward (out of the compressor 
inlet) than rearwards. 

Typical directivity patterns are shown in Fig. 4. 
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1.2 Effect of Walls on Radiation Field 

Since the walls of the wind tunnel are fairly rigid and because 
the acoustic signals we are dealing with are mainly in the audible range, 
the walls of the tunnel behave as good acoustic reflectors* 

The wind tunnel as an acoustic enclosure presents some specific 
peculiarities. In spite of its length, the tunnel cannot be visualized 
from an acoustic point of view as an infinitely long duct to a sound 
source located in the test section of the tunnel. Neither is the tunnel 
an ideal reverberant chamber in the sense that the excited acoustic 
normal modes of the tunnel are in equilibrium with each other. More 
specifically, the combination of sound source plus tunnel does not lead 
to a "diffuse" acoustic field, with a well defined uniform "reverberant" 
energy density. Instead, one expects a large number of modes to be 
excited. Some, whose k vector is along the axis of the tunnel will be 
akin to "travelling, hence nonreverberant modes", others will be standing 
modes with a likelihood of establishing a quasi-reverberant situation. 

We are after all searching for diagnostic techniques that will 
allow us to differentiate between the direct and "reflected" sound in 
the tunnel. The question of paramount importance, therefore, is the 
estimate of a mode "mixing length". A corollary to this estimate is the 
criteria for mixing of lateral or cross modes of the tunnel. 

1.2.1 Quasi-reverberant Regime of Tunne l 

When a simple point source is placed in an enclosure, then under 
steady state conditions, a large number of normal modes are excited. In 
the limit of high frequency and if the source is not near a wall, so that 



15 . 


EL-SUM CONSULTANTS 


NAS 2-6872 


the different standing waves excited have their nodal surfaces distributed 
enough at random that one can use the average value of the intensity of 
these various modes, the reverberant sound intensity in the enclosure 


is 


(15) 


W = 


2 2 
2tt p Q v 
K o o 

2 

a o Y 


( 8 ) 


where Q is the strength of the source 

V is the sound frequency of the source 
y is the average sound absorption of the walls , and 
P o ,a Q are as before the density and velocity of sound, 
respectively. 


Close to the sound source the wave behaves as if the source were 
in free space. Consequently, had one been interested in measuring the 
sound directly radiated from the source, a proximity measurement would 
yield the required information. On the other hand, when the source is 
extended in space proximity measurements are no longer adequate to allow 
one to determine the strength distribution of the sound source. A typical 
example to clarify this point is seen by considering the case of a dipole- 
like source: the nul point of the dipole would not be discernible in 

view of the background reverberant signal. It is clear, therefore, that 
some other method must be designed to map the source distribution. 

But before we actually discuss an appropriate diagnostic technique 
it is important to establish whether it is indeed possible for a rever- 
berant field to exist in an infinitely long conduit (highly idealized 
wind tunnel) . To simplify the mathematical calculations we demonstrate 
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the existence of this reverberant field on a conduit of circular cross- 
section. The conclusions can be extended to a conduit with a different 
cross section, but obviously the criteria are bound to be numerically, 
but not qualitatively, different. 

We begin by calculating the sound pressure field inside a conduit 
with rigid walls and due to a simple source placed on the axis of the 
conduit. It is readily shown from elementary acoustic theory, that when 
the source has unit strength the field is given by 


2 


J (8 r) 

o n 


exp 


~ n 8_ (8a) 

n=U no n 




(9) 


where k is the wave number 2 tt/A 

($ n a) is the n th root of J^(£a) = 0, and 
a is the radius of the conduit. 

The field is described by means of a cylindrical coordinate system whose 

polar axis coincides with the axis of the conduit. Since the source is 

* 

symmetrically located on the axis, there are no normal modes with azimuthal 
dependence. 

Now Eq. (9) indicates that any mode of order n propagates whenever 
k > 3 n> otherwise it is damped. The lowest mode for B q * 0 corresponds 
to the dominant plane wave mode. Suppose we consider all frequencies 
with k < 3^, Consequently, all modes higher than the fundamental are 
attenuated. Figure 5 shows graphically the evolution of the pressure 
distribution from that due to a point source to that of a plane wave. 

Now, in the interim region from the origin to a length l which is 
determined on page 18, the higher order modes do not propagate, and 
classical theory indicates that the energy associated with these modes 




POINT 0.2 0.4 0.6 0.8 1.0 1.2 1.4 

SOURCE 

2/ X 


Fig. 5, Contours of equal acoustic pressure amplitude in a circular conduit. Notice that beyond z/A 
plane wave begins to emerge. 
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represent a reactive (nondissipative) energy. Actually, these modes are 
trapped and bounce back and forth off the walls of the conduit. But this 
is the usual situation for a reverberant field. Consequently, a rever- 
berant field can be set up in a conduit, but it is restricted to the band 
of frequencies defined by k < within the characteristic length Z (a 
function of k) . 

An estimate of the length ,, £ ,t (a distance at which plane wave emerges) 
can be had by considering the exponent of Eq . (9). Manipulation of this 
quantity leads to 




(|3 a) 2 (b 
n a 


Since (g^a) = 3.831, the critical value for a/X is 


,_a, 3.831 

V ■ 2 IT 


0.61 


( 10 ) 


If we require, on the other hand, for the n=l mode to decay by a factor e 
over a wavelength, then we must have 


a = 3.831 

Jl + ( 2 tt ) 2 


.602 


(ID 


The characteristic length Z is then 

1 = (. 6 O 2 ) 


( 12 ) 


i.e., this length is of the order of a diameter. 

The above estimate for the characteris tic length is conserva- 
tive since it was reckoned on the basis of a symmetrically located 
source on the axis. If the source is off axis, normal modes with azi- 
muthal dependence appear in Eq . (9) and the corresponding length l becomes 
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longer. Furthermore, if instead of a simple source we had a dipole 
source, then a larger number of modes are excited. 

The conclusion of these arguments is that in spite of the conduit 
like behavior, it is expected that the sound field in the test section 
of the wind tunnel to behave like a fully developed reverberant field. 

As we saw earlier, the sound source is approximately 8 to 10 turbojet 
nozzle diameters in extension and is equivalent to a mixture of lateral 
and longitudinal quadrupoles . Since the nozzle diameter is much smaller 
than the tunnel characteristic dimensions (40x80 ft) , the sound field in 
the vicinity of the source is definitely reverberant in behavior. This 
does not apply, however, to the whole tunnel. We thus like to refer to 
the field as quasi-reverberant . 

Before we close this phase of the discussion, we will make a few 
remarks concerning the mixing of cross-modes. The exciting of cross 
modes in the wind tunnel is very likely because of the rotating fluctua- 
ting lift forces on the blades of the compressor. The corresponding 
source term becomes a very efficient exciter of the higher order so-called 
"swirling" modes. This question is of importance to this study because 
any mixing of these modes would tend to accelerate the onset of rever- 
beration conditions. We discussed above the situation for a jet alone. 
Now, when the compressor noise is introduced in the discussion we have to 
take into consideration the fact that the large number of coherent 
fluctuations associated with these swirling modes are in fact true 
sources; yet they could be mistaken for random reverberant noise. 

One convenient way of treating this dilemma is to rely on a cor- 
relation for the spectral density of these modes. Bolleter and Crocker 
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have studied this question at length and they were able to show that these 
higher order modes can interfere constructively or destructively depending 
on the frequency and the distance from the source. An axial interference 
"wavelength” L between an "mn" and a "pq" mode (the two numbers, refer to 
azimuthal and radial nodes order) can thus be found and which is given by 

1 ■ r 2 ]i/2 X r C13) 

L 1 - J - L 1 ' < V“ ) J 

here 03 , o> are the resonance frequencies of the inodes and A is the 

mn pq 

wave length in free space. The presence of ,T L ,f is an indication that 
such fluctuations are indeed associated with sources and do not represent 
a spurious noise signal. 
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PART II 

THE DETECTION OF DIRECT SOUND FROM A 
SOURCE IN A REVERBERANT BACKGROUND 

2.1 Introduction 

As pointed out in Part I, the measurement of the direct sound 
radiated from a source in an enclosure can be done by sampling the 
sound field in the immediate vicinity of the source provided that the 
acoustic strength of the source is much higher than the equivalent 
strength of the reverberant sound in the background. If this criterion 
is not fulfilled, a correlation detection technique has to be used to be 
able to differentiate between the direct and background reflected sound. 

Now the use of highly directional microphones allows one to measure 
the sound close to the source and this without having to place the micro- 
phone close to the source. We will show later in this part of the 

report that if a correlation procedure is to be used then the use of 
directional microphones has the very great advantage of enhancing the 
signal-to-noise figure of merit. Directional microphones are thus highly 
desirable instruments in measuring the direct sound. In the following 
sections we first discuss the design of directional microphones and in a 
subsequent section we outline a correlation detection method. 

2.2 Directional Microphones 

An extensive amount of investigation has been done towards develop- 
ing a highly directional sound detector. The basic concept used has been 
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one of preferential wave interference to reinforce the reception in a 
particular direction more than in an another. Two methods have been used 
with success. In the first method, the sound collected from different 
places on a line is chanelled to a single microphone. Several devices 
have been employed^^ , ^ , ^ > ^^ to adjust the acoustic delays in the 
separate channels so that interference between the signals from these 
channels would lead to the required directivity. A typical example of 

this class of microphones is the so-called "tubular or line microphone". 

( 21 ) 

In the second method, several microphones are arranged in an array 
(linear or other) and the electrical output is mixed to allow the phase 
interference necessary to provide the directional response. Experience 
acquired in the design of directional antennas has been extended to the 
construction of equivalent acoustic detectors. We now proceed to discuss 
the design of these two types of microphones. 


2,2.1 Tubular Microphones . 

These consist of a microphone constructed out of a number of 
small tubes with the open end, as pick up points, equally spaced on a 
line. The other end is terminated in a velocity or pressure microphone. 
No matter which termination is used it is imperative that the small 
tubes are terminated into the correct acoustic impedance (pC) to avoid 
the excitation of standing waves in the small tubes. It is readily 
shown that the directivity R Q of such a microphone, Figure 6, is given 
by [Ref. 18] 


R 0 . 


. Sl "[f 


(£ - £ cos 6) 


TT 


j (£ - % cos 0) 


(14) 
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Fig. 6. Line or tubular microphone. Useful directivity on the line 
axis. This microphone consists of a large number of small 
pipes arranged in a line with the distance from the opening 
of each pipe to the common junction decreasing in equal steps. 
The system is terminated in a ribbon element and an acoustical 
resistance. Another alternative is to terminate the system in 
a pressure microphone. No damped pipe is then needed. 


In the above expression Si is the length of the line. The directivity for 
a large n umb er of small pipes (n °°) is independent of n. The corres- 
ponding angular dependence is shown in Fig. 7. 



Fig. 7. The directional characteristics of the microphone shown in 
Fig. 6 as a function of the ratio of the length of the line 
to the wavelength. The polar graph depicts the output, in 
volts, as a function of the angle, in degrees . 
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In order to effectively increase i and hence improve the directivity, 
Olsor/ 18 ^ has suggested to add a delay in each line. The amount of delay 
added is proportional to the distance from the pipe opening to the pickup 
point nearest the common termination (Fig. 8). 



Fig. 8. Tubular microphone with acoustic delays to effectively increase 
the length 


As one would expect it, the same directivity factor is used except that 
the length is now increased. In fact, one obtains 


sm 


R e = 


y a - 1 cos e + d)i 


LA 


TT 


(S, - l cos 0 + d) 


(15) 


Here again Rq has been calculated on the basis of n + <*>, The behavior of 
Eq. (15) is shown in graphical form in Fig. 9, 

The improved directivity of the microphone obviously does not 
represent the optimum that can be obtained with the present-day state of 
the art of antenna design. Furthermore it is seen that the sharpness of 
the directivity is very dependent on the ratio (&/A). Since at the low 
frequencies of the order of 100 Hz \ - 3 meters , it is seen that a line 
microphone would be unusually unwieldy. Consequently a search for another 
diagnostic tool is desirable. 
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Fig. 9. Effective directivity of line microphone with delay line as 
shown in Fig. 8. 

2.2.2 Linear Arrays 

As mentioned above a linear array of identical (directional) micro- 
phones will have a response not unlike that of an ideal line or tubular 
microphone. In both cases, the phase mixing takes place. In one case 
it is with acoustic signals while in the other with the electrical 
(counterpart) signals. The advantage of the linear array , however, is 
that one does not have to be concerned with termination of a line as is 
the situation in the tubular microphone, where poor termination will 
cause mismatch and phase distortions. The drawback of the linear array 
is that one has to use several electro-acoustic transducers, hence they 
have to be compensated to make them identical. Furthermore, the cost is 
higher since one needs M n" transducers. 

The directivity pattern of a linear array can be optimized by 

( 22 ) 

means of a number of methods which have been developed over the 
years (Shading, Dolph-Tchebychef f correction, etc.). All these corrective 
procedures are much more easily adaptable to the electric circuits of the 
transducers, than to the all-acoustic network, so that there is an overall 
benefit in using an array of transducers . 
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Briefly summarized, the directivity of the array can be optimized 

on the basis of equal minor lobes pattern, by adjusting the spacing between 

the transducers, by altering the sensitivity of the transducers or by 

changing both. The other variables one has to consider are the number of 

transducers and the overall length of the array. 

Now, Dolph^ 2 ^ and Pritchard** 2 ^ have shown that the sharpest 

major lobe of a directivity pattern due to a linear array of equally 

spaced point elements is achieved when the elements sensitivity are such 

that all minor lobes in the pattern have the same relative amplitudes. 

It can also be shown that If the major half-width is specified, a maximum 

allowable spacing Is uniquely defined independent of the number of elements 

or the specified minor lobe amplitude. The functional relation between 

(24) 

the half-width, 0 , and the spacing d was given by Pritchard to be 

■L* max 


e L = sin 


, 1 - (d/X) 

-1 max 


(d/A) 


max 


(16) 


If the spacing is made greater than (d/A) x , then the minor lobes no 
longer remain equal, neither is their relative amplitude small by compari- 
son to the major lobe. It is clear, of course, that the sensitivity of 
the various elements on the area has to satisfy a dependence derived by 
means of a Tchebycheff optimizing technique. Indeed, if we use (2M+1) 
transducers equally spaced, so that one transducer is located in the 
center of the array, let us define the normalized sensitivity of the 
individual transducers by b^. The directivity pattern function is 
expressed by 

K ou , , (0) = b + 2b- cos* + ... + 2b w cos M<£ 
zMfl o 1 M 


( 17 ) 
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where <f> = kd sin 0 . The conventional Tchebychef f procedure consists in 
first introducing the variable 


X = COS (f), 

the directivity of Eq. (17) is now rewritten as ^2M+1^ X ^* T ^e re< l u: ^ re ^ 
directivity pattern (i.e., one with equal minor lobes) is now found 
immediately from 

K 2»H (X) " LT M (2) (18) 

where z - (Ax + B) and T__ is the Tchebycheff polynomial of degree M. In 

M 

the above expression L is the relative amplitude of the minor lobes 
(which are all equal) and A, B are constants. These constants are ob- 
tained by means of the following simple procedure. Let z ^ be the value 
of z which makes T.,(z) = l/L, then 

A «* (z + l)/2 and 

Orl 

B = (, oH - l)/2 

The values of A, B and corresponding for different M and L have been 
tabulated by Pritchard (Table I). In fact for L = 1/8 and M = 2 (i.e., 

5 transducers) , we have 

b - 0.281 150 
o 

b 1 = 0.230 032 
b 2 = 0.129 373 

Now if we select 0_ of Eq. (16) to be about 10°, then (d/X) ~ 0.843. 

L max 

Suppose we decide to have the maximum (d/X) = 0.750 < (0.843) then the 

max 

directivity pattern is as shown in Fig. 10. 
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TABLE I 

NUMERICAL VALUES OF SENSITIVITY COEFFICIENTS 
VALID FOR (^ _> j (After Pritchard, Ref. 23) 


L dB 

b o 

b l 

b 2 

-39.7 

0.340 206 

0.247 423 

0.082 474 

-34.5 

0.328 372 

0.245 290 

0.090 524 

-29.8 

0.314 516 

0.241 936 

0.100 807 

-24.6 

0.294 118 

0.235 294 

0.117 647 

-22.0 

0.281 150 

0.230 032 

0.129 393 

-18.8 

0.262 673 

0.221 198 

0.147 465 

-16.9 

0.250 000 

0.214 286 

0.160 714 

-14.8 

0.233 577 

0.204 380 

0,178 832 

-13.8 

0.225 513 

0.199 154 

0.188 090 

-12.8 

0.216 295 

0.192 912 

0.198 940 

-12.0 

0.208 450 

0.187 387 

0.208 388 

-11.7 

0.205 627 

0.185 354 

0.211 833 

-10,6 

0.193 090 

0.176 053 

0.227 402 


Now, the linear array is needed to be used over the entire audio 
frequency range. Since this spectrum covers three octaves, the way we 
propose to build the array is to use one of the following alternatives. 

1. Use 5 microphones whose sensitivities are weighted as indicated 
by Eq. (19). Between 0-100 Hz use a spacing of d = 0.75A = (0.75x3.4) = 
2.54 m. Between 100-1000 Hz use a spacing of .254 m. Between 1000- 
10,000 Hz use a spacing of 0.0254 m. The way to achieve this arrangement 
is to have the microphones mounted on a track and the distances then 
adjusted from outside by a motorized system. 
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The overall extension of the microphone is then 4x2.54 - 10.16 m or 
30 ft. roughly in the maximum position. It is clear that with this arrange- 
ment we make the measurements in three steps, corresponding to three adjustments. 

2. If it is decided to carry all measurements simultaneously, we use 
13 microphones, spaced in three batches. The microphones are then connected 
as shown below in Fig. 11. 



Fig. 11. Connections for 13 microphones arranged in three batches. Each 
batch has 5 microphones for the three frequency ranges 0-100 Hz, 
100-1000 Hz or 1000-10,000 Hz. Microphones are separated by 2.54 m 
in the first batch, 0.254 m in the second and 0.025 m in the third. 

It is to be understood that the weighting functions for the microphones must 
not be frequency dependent . This means that the adjustment is best done with 
purely resistive attenuating circuits. A typical arrangement is then as 
shown in Fig. 12. The criteria as to whether we use scheme (l) or (2) will 
depend on the cost and/or the need to have the measurements done simulta- 
neously or not. This, of course, is set by the cost of running the wind 

tunnel, compared to the cost of microphones and circuitry. 

We cannot complete this discussion without considering another 

method that has been used with considerable success in getting directivity 
through a correlation concept. 
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MICROPHONE 

5 4 12 3 



Fig. 12. Typical block circuit diagram for one of the batches of Fig. 11. 
2 . 3 Correlation Detectors 

The use of a correlator with a two-element array of directional 

receivers is one of the simplest arrangements of correlation detectors 

(25) 

which has been used in radio astronomy investigations. One of the 

chief advantages of such a system is that it is very insensitive to 
background noise from a distributed source (such as reverberant noise) 
allowing one to use greater recorder sensitivity than in the conventional 
sum-and-square law detector systems, thus making this technique specially 
suited to the detection of weak signals* 

The output of the detector or correlator will, in general, consist 
of several components: a steady or d.c. component due to the signal; a 

steady component due to background noise; a fluctuating component due to 
the signal and a fluctuating component due to background noise. If the 
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third contribution is small compared to the last (i.e., small signal to 
noise input to detectors) we can neglect it. If we use a directional 
receiving system then in principle the signal (here the direct sound) 
•should fall to zero as the signal source is removed from the principal 
axis of the receiving system. The d.c. output from the receiver thus 
corresponds to the background noise and this part could be biased off. 

The effective response pattern of the system is improved. Indeed 
if the receiving pattern of one microphone is K(0) , the relative response 
of the correlation system is K (6) p(d/c sin0) where p(r) is the normal- 
ized autocorrelation function of the signal. The delay T * d/c sin0. 

If one had used two tubular microphones, then the system can 
indeed pick up the direct sound and discriminate against the reverberant 
noise. Inspection of Figs. 7 and 9 indicate that the directivity of 
this microphone is poor by comparison to that of the linear array. 

2.4 Concluding Remarks on Direct Detection of Sound Source in a 

Reverberant Background 

The linear array appears to be the simplest of the various schemes 
we have discussed. It also displays the sharpest beam. For this reason 
we have decided to use it for measuring sound directly radiated from 
source. 

To be sure, more complex (two dimensional) arrays can be found 
which would display an improved characteristic. We propose first to start 
with the linear array to demonstrate the feasibility of the technique. 

So far we did not discuss the manner with which the array 

will be steered. In view of the physical size, it is clear that it is 
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much easier to steer it electronically. The procedure is quite straight- 
forward. 

Finally, as we pointed out above, two alternatives were suggested: 

(a) variable array length, and (b) fixed length many microphones used. 

We feel that the second alternative may be more costly in adding to the 
numbers of electro-acoustic channels, but the increased speed of perform- 
ing the experiment more than outweighs the cost of the added equipment, 
through saving of personnel time and cost of operating the tunnel for 
longer periods. 

The design that ARC approved is that for a crossed array, five 
elements per array. This is described in detail in part IV and Appendix 


B of this report. 
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PART III 

OPTICAL DIAGNOSTIC METHODS 


3.1 Introduction 

The optical methods for the diagnosis of noise from a jet flow, 
which we proceed to expound, are based on the use of coherent optics 
and opto-acoustical techniques. These methods enable us to map the 
fluctuating (turbulent or coherent) part of the flow and hence measure 
small density variations and check the theories of mechanical to 
acoustical energy conversion. Furthermore, localization, spectrum 
analysis, background discrimination and optical processing of a wide- 
band signal will be facilitated. 

These methods are grouped under holographic interferometry, 
acoustical holography and optical processing of wideband acoustic 
signals . 

3,. 2 Holographic Interferometry 

Fundamentals of holography are well known. ' It is based on 

recording the interference between a known and an unknown wavefront. The 
latter can be reconstructed when the original recording (the- hologram) 
is reilluminated by the known reference wave. 

Holographic interferometry depends upon the interference between 
two reconstructed wavefronts. These wavefronts may have been recorded 
on the same hologram (double exposure) or on two different holograms, 

( 28 ) 

reconstructed in such a manner that the reconstructed waves coincide. 
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Slight deviation between the two reconstructed waves is revealed in the 
form of three-dimensional interference fringes superimposed on the re- 
constructed image. The fringe data can be reduced to density informat ion * 

(29) (30) 

■as was done by Collins v who made use of Maldonado^ ' line integral 
inversion scheme for the development of plasma emissivity . In this 
scheme, the fringe number functions that describe the interferometric 
data g(y T ,£) is used to obtain the density-related function f(x,y) by 
inverting the integral 

x ? 

2 

g(y\£) = Q f f(x»y) dx’ (20) 

x T 

l 

where 

Q = p„ Mp s (2D 

f(x,y) = p(x,y)/p oj - 1 (22) 

n = 1 + gp/p s (23) 

p g = reference density at p°C and 760 mm Hg 

8 - dimensionless constant related to Glads tone-Dale constant 

n = index of refraction 


Collins and his co-workers reported that the fringe data may then 

be reduced to density information by expressing the density in a series 

of orthogonal polynomials and inverting the fringe-density integral 

equation by a computer program in the manner described in Ref. 29. Their 

scheme in brief consisted of using Q-switched laser to take double ex- 

-9 

posure holograms (20x10 sec. exposure) to freeze the jet flow. The two 
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exposures were taken before and after the jet was turned on. The fringes 
in the reconstructed images were then analyzed and the computer program 
was used to determine the density distribution. 

(31) 

A different arrangement was used by Jahodo to map the auto- 
correlation of the density. The arrangement is based on taking an in- 
focus-diffused hologram of the flow. The density distribution is revealed 
as light intensity variation within the field of view. 

3.3 Acoustical Holography 

3.3.1 Direct Acoustical Holography 

The acoustical field may be mapped into an optical one by making 
an acoustical hologram which is reconstructed optically. There are 
several ways of making such a hologram, which are summarized in Ref. 26. 

In this project, the information obtained by using the linear acoustical 
detectors (Section 2) may be used to construct a hologram of any region 
of the jet flow. This hologram, reconstructed optically, will then 
reveal the coherent sound sources as visible images in the right location 
with the correct intensity. In brief, this approach is rather another 
way to analyze the data obtained by the linear detectors array described 
in Section 2. 

3.3.2 Instantaneous Acoustical Holography 

Two superimposed pulsed holograms of the jet flow are made within 
a quarter of an acoustical cycle. The reconstruction of such a compound 
hologram reconstructs the acoustic field optically, indicating both the 
location and intensity of the coherent acoustic sources within the field 
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of view. This is an all optical technique where no acoustical trans- 
ducers are needed. However, it is limited to the detection of purely 
coherent sources of known frequency; furthermore, it is rather costly, 

•due to the need for strong double pulsed laser and good optical components 
capable of spreading the laser beam to a large extent (100 cm or more). 

3.3.3 Spectral Analysis and Source Localization From Acoustical 

Holograms 

When the acoustic field is mixed with a reference signal of fixed 
frequency f, and scanned by an acoustical detector with a constant 
velocity v, the resultant hologram reconstructs the acoustic field at an 
offset angle 0 given by 

Sin 9 = j • (24) 

where C is the sound velocity and Af is the frequency difference between 
the unknown source and the reference signal. Thus, it is easy to see 
that discrete sources of different frequencies will be reconstructed 
at different angles 0. In other words with this simple approach noise 
sources may be localized and their fundamental frequencies determined. 

Each technique, however, is valuable for stationary sources having 
single characteristic frequencies. Information on other types of sources, 

( 0 fa} 

however, can be obtained but with difficulty and less accuracy. A 

more serious drawback, though, is the needed stability of the acoustic 
field during the recording of the whole field. Again, an acoustical 
detector is needed and the linear array described in Section 2 can be 
used to form such a hologram. Better still, two crossed linear arrays, 
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if available, will be very useful and time saving. 

3.4 Optical Processing of Wideband Acoustic Signals 

An elegant method to determine the power spectrum of acoustic 

(32) 

signals, using coherent optical techniques, was developed by Markevitch. 

The principle of the technique is simply to transfer the recorded acoustic 
signal to a two dimensional photographic transparency in the form of an 
appropriate raster. An optical Fourier transform of the transparency 
will then reveal the spectrum of the recorded signal. Spectral resolu- 
tion of 0.5 Hz is easily obtained and 0.05 Hz was demonstrated. Signals 
of 30 dB below the background noise was detected and 50 dB below back- 
ground noise is hoped for. 

The equipment to carry out the analysis is rather very expensive 
to build, especially just for this project. Instead, already existing 
equipment can be used to analyze a tape recording of the noise by special 
arrangement . 

3. 5 Concluding Remarks on Optical Diagnostic Methods 

In addition to the linear detectors array described in Part II, 
availability of (a) pulsed laser, (b) CW laser and (c) appropriate optical 
components and mounts are needed to make and reconstruct either optical 
or acoustical holograms. The analysis of the data may be done directly 
by optical means and by computer. 

By special arrangement, the detected signal from the linear array 
may be tape recorded and analyzed by the coherent optical processes. 

Measurement of the density fluctuations in the flow by other 
optical means are discussed in Secs. 3.2.3 to 3.2.6 of Appendix A. 
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PART IV 

DESIGN CHARACTERISTIC OF A HIGH DIRECTIVITY 
ACOUSTIC ANTENNA 

4.1 Basic Considerations 

The response of a planar microphone whose transduction efficiency is 
w(x, y) can be shown to be proportional to 

Jf w(x, y)exp ik(x cos a + y cos 6)dxdy (25) 


where the interpretation of the various geometric symbols are readily found 
from Fig. 13. 




Incident Wave 



Fig. 13 Geometry of 

Antenna Array. 


For a distribution of point microphones, the response is proportional to a 
quantity which reduces the double integral to a double summation, viz.. 


V / w w exp ik(x cos a + y cos (3) 
L L n m « « 


n m 


n 


n 


(26) 


where w , w are the individual sensitivities of the microphones and n, m 
n m 


are the number of microphones in the x and y directions respectively. 
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4.2 Steering of Antenna 

When the antenna is tilted by a (or 6 ) there is an acoustic phase 
2xfd 

difference of — cos a (for 3=0, say) between neighbouring micro- 

c 

phones, where d^ = distance between microphones in the x-direction, f is 
the frequency and c the velocity of sound. 

In order to steer the antenna by electronic means , one has to introduce a 

relative phase shift between the outputs of the different microphones. To 

produce an effective rotation of the (x-arra) by an angle a , the relative 

2TTfd 

X 

electronic phase shift must be ( cos a) - Similarly a rotation of the 

c 2TTfd 

(y-arm) by an angle S would require a phase shift of — cos $ . 

The important design parameters which have to be decided upon in order to 
provide the design criteria for the electronic circuitry are: the number of 

microphones to be used and the separation distance d between the microphones. 

Of course the separation in the x direction could be different from that in 
the y-direction* For equal sharpness in the x and y direction the separations 
d^ = d^ = d say. The actual magnitude of d will depend on the Tchebycheff 
design criteria which was discussed in Sec- 2-2-2. 

The number of microphones to be used is set by the directivity sharpness 
desired together with cost considerations. Figure 14 shows 

that the gain in sharpness from a (15°) to a (5°) sharp beam would entail an 
increase in the number of microphones from 5 to 15 per direction (i.e. from 9 
to 29 total for the same discrimination of 26 db.) This gain in directivity 
represents an appreciable increase in capital investment which is not warranted., 
Since the total number of microphones available was thirteen which would have 



MINOR LOBE AMPLITUDE RELATIVE TO UNITY, DB 



Fig* 14 Minor-Lobe Amplitude as a Function of the Length of the Array. The angle 0 refers 
to the half angle of the beam. The actual number of microphones can be deduced from 
the ratio d/A . 
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yielded an array of 7 microphones per direction, it was decided to use a total 
of 9 microphones (i.e. 5 per direction). This decision was justified on the 
grounds that the gain from 25 dB to 36 dB microphone was not appreciable. 


As was discussed earlier, Section 2.2, the Tchebycheff correction is 

d/A dependent. Concentrating on the frequency range 100 Hz to 1 kHz with a 
d/A = 0.75, the maximum electronic phase shift will occur at a frequency of 
1 kHz. Denoting the various microphones by a_^, a^, etc. — and b _ 2 ' 
etc. ..., we have typically (d = 25.4 cms) 

TABLE II 


Station Phase Shift 


Phase Shift for a = 0° 
at f = 1 kHz 


-2 


-1 


4.7 f 

10 ' 


cos a 


0 


2 ( 


+1 


3 ( 


+2 


4 ( 


0 

4.7 radians 

2 x 4.7 

3 x 4.7 

4 x 4.7 


The phase shifters are all-pass networks designed to provide a phase 
shift of 4.7 radians at 1 kHz. A typical such network is shown in Fig. 15. 

The phase shift for an all-pass network is a linear function of the frequency 
and is in fact proportional to the time constant RC for the circuit. But 
the phase per unit (as shown in Fig. 13) is limited to about 40°. Consequently 
8 similar units connected in cascade are needed to provide the required phase 
shift, and one can express the electric phase shift 6 by 


0 = lu F (RC) 


( 27 ) 


where F is a complicated function of RC and o> = 2nf . 
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Now, we equate the acoustic phase shift to the electric phase shift, i.e., 

, , 2ufd 

u)F(RC) - — — cos a 

or 

cos a = cos *[F(RC) ^ ] (28) 

We thus notice that the steering angle a is a function of (RC) the electric 
time constant of the circuit and of the separation between the microphones, but 
is independent of the frequency. 

Consequently, the steering is obtained by changing the value of the resistance 
R , which in the circuit of Fig. 15 corresponds to the resistance of the FET . 

By using a simple feedback control circuit, all the FET's in one phase shifting 
block can be made to change together by changing the voltage on the FET . There 
will be, of course, individual variations for the phase shift produced by the 
different units because of the unavoidable variations of the characteristics 
Of the FET . By trial and error of hundreds of FET , the (2x8x9) = 144 FET's 
with the best match were selected. 

The calibration of the instrument i.e., dial reading vs. average electric 
phase shift is shown in Fig. 16, while the transformation of Fig. 16 by means of 
Eq.(28) gives the steering angle for different values of d (compared to 
j 3 

(2- = —) at 500 cps . ) Fig. 17. 

Because the FET's cannot be made to reach an infinite resistance (by removing 
the applied voltage), it is not possible to satisfy F(RC) -*• 0 , i.e., the zero 
phase shift point cannot be attained with the FET's alone. Consequently, relays 
have been added to by-pass the outputs from the microphones to be added via the 
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operational amplifiers without any phase. The arrangement shown in Fig. 18 
is self-explanatory. 


4.3 Beam Sharpening 

The shading of the microphone sensitivities has been described in details 
in Section 2.2.2. For 5 transducers per arm, the weight of the coefficients 

are as indicated below 

TABLE III 


b „ 

= b ^ 

= 0.46 = a „ = 

a 

-2 

+2 

-2 

+2 

b , 

= b , 

= 0.82 = a = 

a , 

-1 

+1 


+1 

b 

o 

- a 

o 

= 1.0 



The Tchebycheff corrections reduces the minor lobes but does not sharpen the 
main lobe. 


A simple procedure to achieve this sharpening is to use an elementary 

(integrated circuit) multiplier. The multiplier which effectively performs 

the operation (£ a ) x (£ b ) sharpens the beam by 6 db at the half power 
n m 

point. The disadvantage of the multiplier is that its output is at double the 
frequency. It is possible, however, to use an automatic gain control amplifier 
to restore the output at the same frequency, but with an amplitude proportional 
to ( £ a^) x (£ b^) thus preserving the beam sharpness. The multiplication 
scheme obviously necessitates tuning the outputs of each array to a single 
frequency . 




Pig, 18 Interconnection for M A n Channel. An identical circuit is used for the B 
channel. Notice the circuit used for the adding operational amplifiers so as 
to yield unity gain. 
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4.4 Testing of Antenna 

A series of tests were performed, before delivery of th.e antenna control, 
to check on the operation of the system. These are described below. 

1. Phase Shift . The individual phase shifts produced by the 8 delay 
lines (phase shifters) have been checked. Typical phase shifts 
produced have been shown in Fig. 19. 

2. Linearity of Delay Lines . Keeping any one delay line fixed at a 
particular frequency, the amplitude of the input to the delay line 
was varied until a 1% occurred at the output. This limiting 
amplitude was found to be 50 mV. It is important to stress that 
all operational amplifiers gains have been adjusted so that their 
gain was always unity. The amplitude of the signal at the micro- 
phone plugs is therefore identical to that at the delay line input 
even though it goes through a series of amplifiers. 

3. Linearity of Amplifiers . By connecting all the test points in 

parallel, one can effectively add the inputs of the amplifiers. 

Using the "A 11 channels alone and "B" channels separately it was 

checked that the signal J a (or £ b ) did in fact increase 

n 

with n , when the gain of all the channels as fixed by the potentio- 
meters was made identical. 

4. Linearity of Multipliers . Setting the phase shifts of the delay 

lines to zero, it was checked that the output () a ) x (j b ) 

u n u n 

increased linearly with either ([ a^) or (£ b^) . Furthermore, 

it was checked that the product went to zero if either (£ a^) or 

(J b ) went to zero, 
n 
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Fig. 19 Scatter in the characteristics of the phase shifters for the 

"A" channel. The curve drawn is the average for all the delay 
lines. 
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It was also found that the multipliers begin to distort when the 
signal length at the input of the system exceeds 30 mV. For this 
reason the signal from the microphone should be kept to leas than 30 mV, 

4.5 Evaluation of the Antenna and its Control 
4.5.1 Procedure 

With the cooperation of Mr. Paul Soderman of Ames Research Center, 
the following tests were made at ARC anechoic chamber to evaluate the per- 
formance of the acoustical antenna and its electronic steering control. 

Nine omnidirectional microphones were mounted on two crossed arms 
A & B, 5 microphones per arm (Fig. 13), 45.7 cm (18") apart. This antenna 
was placed at one end of the anechoic chamber, facing a loudspeaker at the 
other end, about 6 m (20 ft) away. The antenna could be rotated around its 
vertical axis through 360°, in synchronism with a chart recorder. The con- 
nections and shading of the microphones followed the instruction manual, 
reproduced in Appendix B. 

(a) In one series of tests, an acoustic frequency was selected; no 
shading was made; the controls of the electronic steering of A and B were 
switched off to simulate a (0,0) position of the antenna; the antenna was 
then rotated mechanically in synchronism with the chart recorder which, plots 
the output. This was repeated at various other selected frequencies. 

(b) The same set of runs were made after the microphones have been 
shaded. 

(c) A third set of runs were made (with and without shading) when the 
distance between the microphones was changed to 22.9 cm (9"). 

(d) Another test was made which confirmed that the output of each arm 
alone is the sum of the outputs of the microphones on such arm. 
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(e) To confirm the insensitivity of the equipment to frequency, an 
electronic signal, simulating the acoustical one, was applied to all the 
microphones. The output then was seen to be the same regardless of the 
frequency of the input signal, or the steering controls. 

(f) Finally a few runs were made to test the electronic steering, by 
leaving the antenna in the stationary position and plotting point by point 
the output versus the steering angle as read from the position of dial A 
alone. (Changing dial A alone simulates rotation of the antenna around 

a vertical axis.) 

4.5.2 Results 

Fifteen samples of the results of tests (a) to (d) of the previous 
section, Sec. 4.5.1, are enclosed. A summary of the analysis of these runs 
is tabulated below (Table IV). From these results, three main observations 
are made: 

1. The antenna is indeed directional. 

2. Shading improves the discrimination between the main and secondary 
lobes, but also widens the angular field of view. 

3. When the antenna is adjusted for an optimum d/X at a certain 
frequency, lower frequencies can be tolerated more than high frequencies which 
confirms our theoretical finding on page 26. At high frequencies, the amplitudes 
of secondary lobes increase and confuse their discrimination from the main lobe. 
Optimizing ( d/X ) at 500 Hz gives reliable measurements in the range from about 
200 to 750 Hz (see Fig. 22); if (d/X) is optimum at 1000 Hz, measurements will 

be reliable from about 550 Hz to 1400 Hz (see Fig. 22). 

4. Experimental results agree with the theoretical results as evident 
from comparing, for example, angle 0 and the positions and amplitudes of 
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maxima and minima of run it 7 (Fig. 20) with the theoretical calculations 
(Fig. 21). 


TABLE IV 

SUMMARY OF RESULTS OF ANTENNA EVALUATION 


Symbol 

No. 

Frequency 
_f Hz 

Mic. 

Separation 

d_ 

cm (in') 

Shading 

Main 

Angle 

Lobe 
0° at 

Diff. between 
main lobe and 
secondary lobe. 
Amplitudes DB 

Half Width 

6 db drop 

(4) 

560 

4c (18) 

No 

8 

7 

9 

• 

300 

4C (18) 

No 

17 

17 

12 

(1) 

200 

4C (18) 

No 

17 


5 

(1) 


(18) 

No 

8 


5 


560 

4(. (18) 

Yes 

18 

11 

19 

(2) 


4(, (18) 

Yes 

10 

6 

16 

(3) 


4C d 1 2 3 4 * * * 8 ) 

Yes 

36 

35 

19 

(4) 


23 

No 

12 

10 

14 

(3) 

350 

23 < 9 l 

No 

37 

32 

22 

(2) 


23 (^ 

No 

7 

7 

14 



23 ( 9 ) 

Yes 

12 

10 

17 


560 

23 ( 9 ) 

Yes 

30 

21 

17 


■ 

23 ( 9> 

Yes 




(2) 

■ 

i 3 (9) 

Yes 

11 

8 

18 

(1)&(2) 



Yes 

5 

4 

17 

J 


Symbol: 


(1) Measurements are approximate. 

(2) Amplitude of high order lobes reaches that of main lobe. 

(3) High order lobes are indistinguishable. 

(4) Optimum (d/A) = 0.75. 


The results of the electronic steering test (f) described in Sec. 4.5.1 

were not quite conclusive. There was indeed an indication of steering, but 

the high sensitivity of the steering dials A & B at low angles made the 

* 

evaluation rather difficult. Further lengthy tests requiring fixing the dials 


A and/or B at a specific angle and repeating tests (a) to (d) were not done. 


In a new model, the phase shifters are so designed that the amount of cascading 
of the control resistance is a function of the steering angle. This equalizes 
the sensitivities of the dials thoughout the steering range. 
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Fig. 20 Fifteen charts of runs made to evaluate the acoustical antenna 
performance. The main parameters (separation d between the 
microphones, the frequency of the acoustical wave, and whether 
the microphones are shaded or not) are indicated on every 
chart. All the runs are made with the electronic steering 
control dials A&B in the off position simulating a (0,0) 
directivity position. (See Appendix B.) 
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Fig. 21 Acoustical antenna response at 500 Hz (five shaded microphones per 
arm, placed 18" %part) . The output is the product of the two arms 
of the antenna. Solid curve is theoretical; dotted curve is experi- 
mental. Both curves agree for 0 at half maximum (6 db) , and for the 
difference (19 db) between the main and minor lobes. Shift in 
position of high order lobes in the experimental curve and lack of 
cancellation at the minima indicate possible poor phase matching 
between microphones, misalignment of the antenna (not working exactly 
in the far field) , possible difference between individual microphone 
responses, and existence of white noise. 
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Fig. 22 Drop between the main lobe and secondary lobes as a function of 

frequency for acoustical antenna with microphones 9" *(solid curve) 
and 18" v dotted curve) apart; d/A = 0.75. Also, the semi-angular 
resolution 0° at 6 dB for the two cases are shown. For a minimum 
of 6 dB drop and to avoid confusion with antenna characters, the 
9" microphones separation is operative within 550-1400 Hz, and the 
18" microphones separation is operative within 200-750 Hz. 
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PART V 

CONCLUSION AND RECOMMENDATIONS 

5. 1 Conclusion 

The feasibility of measuring aerodynamic noise in a reverberant environ- 
ment (e.g. , the 40x80 foot wind tunnel at ARC) without disturbing the flow has 
been shown. Of the various techniques which can be used for this measurement, 
purely acoustical methods (specifically, acoustical antennas), were adopted 
because of the availability of many of the needed components. Equipment was 
built and evaluated to carry out such measurements. Other techniques including 
purely optical or opto-acoustical approaches requiring either more sophisticated, 
less available, or more expensive equipment, are still needed to supplement the 
acoustic antenna results for complete mapping of the acoustic field and correlate 
it with the theories of noise generating mechanisms in an aerodynamic flow. In 
brief, the acoustical antenna scheme which has been built is able to: 

(a) localize the noise sources with an angular resolution of about 10° at 6 db, 

(b) determine the frequency and strength of such sources, 

(c) discriminate between the source and reverberant background, 

(d) do all the measurements without disturbing the jet flow, 

(e) steer the acoustic antenna electronically without disturbing the 

spatial position of the antenna, its components, or readjustment of the electronics 

* 

of the control circuit* 

(f) measure various frequency ranges with the same electronic equipment* 

by merely changing the separation between the antenna microphones, or using other 
sets of microphones which are prearranged (Sec, 4.5.2 (3), p. 52 and Sec. 2.2.2 (2), 
p. 30), and 

^Microphones are assumed to have flat characteristics from 0 to +90° only; other- 
wise the tunnel wall behind the antenna should be covered with absorbing material 
to prevent back reflection. 
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(g) do all the measurements without readjusting the electronics because 
of any change in the radiation frequency under study. 

We may add that microphone shading resulted in better directivity response by 
equalizing the minor lobes (Sec. 4.5.2 (2)). The shading coefficients are 
constant as long as the ratio d/X is unchanged. 

Now, we proceed to expound, briefly, our recommendation for future studies 
on this project. 

5.2. Recommendation for Further Studies 

1. A new model of an acoustical antenna, the design of which is basically 

the same as the one that has already been built, is recommended. This new model 
will feature: (a) easier operation; (b) maximizing the performance in the 

sense of obtaining higher discrimination (higher angular resolution), same out- 
put at input frequency, sharpening the main lobe, suppressing minor lobes, 
automatic shading, etc. , using the same number of acoustic transducers (micro- 
phones) if possible; (c) a self-contained portable package which can be used 
indoors or outdoors; (d) near field measurement capability, as well as far field 
measurements; and (e) automatic as well as manual steering. 

2. Measurement should be done to correlate the acoustical measurements 
with the density fluctuations in the flow. Optics techniques, probing the flow 
(looking inside it) are available. These techniques are of two types: 

(a) point by point probing, and 

(b) interferometric and holographic mapping of the acoustic field. Part 
III and Appendix A, Sec. IV, deal with these methods extensively. 

3. Design a new acoustical antenna, having the general capability of 
the already designed one (Model 909B) ; and, in addition, producing an acoustical 
hologram to be analyzed later on, either by computer or optically, to confirm 
the noise identification as well as the density fluctuation, simultaneously. 
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APPENDIX A 

EARLIER STUDY OF THE PROJECT 


The following 29 pages are reproduction of an earlier analysis 
of the project . This analysis were made by the authors and 
their associates. 
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I. INTRODUCTION 
1.0 GENERAL REMARKS 

The basic formalism for the aerodynamic generation of noise in 

( 1 ) 

its broad outlines was developed a few years ago by Lighthill. 

He was able to show that the general equations of fluid 4Y narn i cs ' 

by means of a quasi-linear approach lead to a formal expression 

for the sound generation, wherein the fluctuating Reynolds stresses 

( 2 ) 

play the equivalent role of acoustic sources. Lighthill himself 
and several other investigators ^ ^ ^ after him pointed out that 

the strength of his effective sound sources is in fact the first 
term of a series expansion in the Mach number for the flow. Since 
this series diverges for Mach numbers greater then unity, Lighthill 1 s 
earlier formulation is applicable solely to subsonic flows . 

The mechanism of noise aerodynamically generated by a fluctuating 
supersonic flow is expected to be quite different from the subsonic 
regime. This is not surprising at all, since fluid elements moving 
coherently and with supersonic velocity will radiate collectively 
in a manner reminiscent of the familiar Cerenkov effect. Accordingly, 
a number of semi-heuristic models have been proposed in 

which the acoustic effect of the fluctuating flow is represented 
by virtual sources made up of discrete "eddies’ 1 radiating as moving 
sources (Eddy Mach waves) . 

In light of the brief discussion mentioned above, we will survey 
the details of noise generation in the subsonic case separately 
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from that in the supersonic case. It will be seen that Lighthill s 
theory is a very powerful tool which allows the treatment of any 
aerodynamical ly noise generation - mechanism associated with a 
subsonic flow. 

1.1 THE SUBSONIC REGIME 

We will find it very convenient from a diagnostic point of view 
to discuss two subdivisions of noise generated in this flow regime: 
a) The broad band noise (turbulence) , b) the coherent or narrow 
frequency band. 

1.1.1 noise from turbulent flows 

Since reference to Lighthill' s theory will be made quite often in 
this section, we will summarize some of its salient features. It 
is readily shown (Ref. (1)), that the sound field is described by 


^ P 2 2 

<r at 


a 2 P _ i !m 


_ 2 8x. Sx. 
C o 1 J 


( 1 ) 


where P is the fluctuacting acoustic density, c Q is the velocity 
of sound in the fluid region which is at rest (i.e., outside the 
flow) and T ± ^ stands for (p v i v j + P ij " c 0 ^ 2 P <5^). The 

beauty of Lighthill' s formulation is that it not only takes into 
account the effect of the fluctuating Reynolds’ stress p v ± y , but 
also nonlinearities in pressure and densities. In most cases of 
interest here, the Reynolds' stress are the dominent source terms. 
Consequently, as was argued by Lighthill the right hand side of 
the equation are equivalent to quadrupoles . 
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Now for randomely fluctuating field, the intensity is the 
essential sound field variable. Hence it is found that the inten- 
sity I is derived from 


Kx) 


16 f p 



o 


fr ( vV ( vV ( v 2 k> <y *t> 
JJ |*-x | 2 |*-z | 2 





The domain of integration is carried over the whole region of 

2 2 

the flow. From eq. (2), it is clear that 2— t. . T, 0 

9t 1J at 2 

is akin to a correlation. Hence a contribution will only come 
from the region of the fluid which is moving coherently. If 
one defines a characteristic dimension L within which the velo- 
cities are correlated, one can divide the field into an imaginary 

3 

set of non-overlapping regions of volume of the order L . The 
intensity given by (2) is now visualized as the contribution of 
a distribution of point quadrupoles localized in these non- 
overlapping regions. The problem of exactly estimating the strength 
of these quadrupoles requires an a. priori knowledge of the fluc- 
tuating flow field {i.e. turbulent field) . 

By taking the Fourier transform of (2) one obtains an expression 

for the spectral intensity evaluated from a knowledge of the 

3 2 

spectrum of the turbulence (i.e. of — T ). This has 

at 2 

been calculated for a turbulent field with no mean flow (Proud- 
man^, Mawardi^ 8 -). When the flow has a mean flow, the eddies 
the radiation field causes the sound pattern to become 


move, so 
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distorted {from that of a stationary quadrupole) . But it is clear 

that the tensor T. . is strongly affected by the presence of shear. 

X J 

Now the shear in the flow may be brought about either by a layer 
of fluid moving against a fluid at rest (case of a jet) or by a 
layer of fluid moving in the vicinity of a solid boundary. These 
two cases will now be looked into in some detailes 
a) Free Jet 

It was argued by Lighthill that the boundary layer in the jet 
tends to polarize the quadrupoles , thus establishing a stronger 
coherence and a consequent reinforcement of the sound radiated. 

By means of dimensional argument, one finds that the acoustic 
power radiated at a distance is given by 



o 


where the coefficient K is a very slowly varying function of the 
Reynolds number. Typical experimental findings are: 

i) The efficiency of conversion of mechanical to acoustical 

energy follows the predicted value M 5 fairly well, (where M=V/c q ) . 

-4 

Typical values of K are 0.3 to 1 x 10 

ii) The sound spectrum is very broad and has an ill defined 
maximum which seems to vary as V/d. 

iii) The polar distribution of the sound has a "butterfly" type 
configuration in which almost all the sound radiated is encompassed 


in an acute angle with the jet. The theoretical directivity of 
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of (1 - M cos 0) 6 seems to be reasonably verified. Indeed 

Hubbard and Lassiter (9) find a 5:5 power dependence (instead 

(f cm) <7.6 e*> 

of 6) for a 3/4 inch’ A jet, 7.0 for a 3 inch^jet and 6.3 for a 12 

inch* j elf. Westley and Lilly (10) confirm the 6.3 law over a wide 
A 

range of diameters. 

It is interesting to note that after a considerable accumulation 
of experimental observation ^^over the intervening past twenty 
years, there seems little doubt as to the validity of the Light- 
hill's theory for the noise from a subsonic jet. 

b) Effect of Solid Boundary on Noise Production 

Lighthill’s expression ^ 2 * can still be applied here, pro- 
vided that one uses for T^j the fluctuating shear stresses exist- 
ing in the boundary layer in the vicinity of the solid body. Now, 
in a boundary layer, it is the motion in the region of fully 
developed turbulence outside the viscous sublayer which sets a 
distribution of velocity gradients responsible for the shear stress 
at the walls (12) . In a fully turbulent region there is a whole 
assembly of eddies which move with the different convection veloci- 
ties. Since the motion in the region outside the viscous sublayer 
is governed mainly by inertia forces , the frequency scale is- de 

fined by 



( 4 ) 
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V is the mean flow velocity. Phillips was able to show that the 
portion of maximum influence is 1/10 of the thickness of the layer 
and that the acoustic power radiated is approximately 


p 0, 



p o V 


(5) 


K is a slowly varying function of the Reynolds number of the order 

of 10 ^ and M is a mean Mach number, 
c 

It is worthwhile to note that the sixth instead of the eighth 
power dependence for the velocity in (5) is indicative of noise 
radiated from a distribution of dipoles instead of quadrupoles . 
This is not surprizing since a quadrupole in front of a solid 
boundary radiates like a dipole (i.e. the quadrupole plus its 
mirror image) . 


1.1.2 THE COHERENT NOISE CASE 

The generation of acrodynamically produced sound of single fre- 
quency (or of very narrow band) is invariably found wherever some 
solid boundary allows the possiblity of some feedback which tends 
to stabilize the fluctuation in the flow around a particular mode 
of oscillation of the system (edge tone, vortex shedding, etc.). 
There again we find Ligh thill's equation (2) of use, provided we 
use in addition an expression for the acoustic feedback on the flow. 
The essential point of interest from a diagnostic point of view is 
the complete description of the flow mechanism if one can deter- 
mine the spatial distribution of the velocity fluctuations. 
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1.2 THE SUPERSONIC CASE 

The very nature of the acoustic radiation from an element of 
fluid moving with supersonic velocities makes it very unlikely 
to obtain sound generation with narrow band spectra (or coherent) . 
Accordingly only the free jet and the effect of solid boundaries 
will be considered. 

1.2.1 THE FREE JET 

The theory of supersonic jets has been considered theoretically by 
Ffowes Williams (5) and by Ribner (6) . Ffowes Williams has shown 
that for this case the sound emission is analagous to the Mach 
wave emission from supersonically moving eddies, the latter behaving 
very much like solid bodies . The sound source is hence equivalent 
to simple source and the source strength is estimated from a density 
correlation by an expression of the form 



u is of course the mean value of the flow. 

The effect of the shear appears as an amplifying factor. 

(14) 

Now it was pointed out by Nagamatou et. al., that an im- 

portant property of a supersonic jet is that the flow contains 
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a region that is supersonic, a transition region and finally 
a subsonic fully developed turbulent region « 

The noise generation accordingly is different for the different 

regions . In the supersonic region the sound power radiated 

follows a M^ 2 dependence (Ffowes Williams) ; Nagamatou on the 

7 — ot 

other hand, suggests that the formula be replaced by M 
near the jet exit and M 7 “ 6 at the end of the supersonic region. 
Experimental values yield a 6. 2 > $ % 2.4 . It 

appears, however, that the a, 3 coefficients are sensitive 
functions of temperature, velocities and atomic weights of the 
gas so that the theory is still far from being settled. 

In summary we find the table shown below very helpful for the 


discussion. 


AUTHOR 

FLOW REGIME 

NATURE OF SOURCE 

IMPORTANT 

PARAMETER 

Lighthill 

subsonic jet 

quadrupole 

velocity 

correlation 

Phillips 

subsonic 
boundary layer 

dipole 

velocity and/or 
correlation 

Ffowes 
Williams , 
Ribner 
Nagamatou 
et. al, 

supersonic jet 

monopole 

density 

correlation 


supersonic 
boundary layer 

i — — 

dipoles 

J 

pressure 

correlation 
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II. FREE FIELD VS. REVERBERANT FIELD MEASUREMENTS 
The radiation field from point sources (monopoles, dipoles, etc.) 
in free space is well known and yields well known directivity 
pattern as well as divergence dependence ( 1/r 01 ; a =1 for 
simple source, a =3 for dipole, etc.) . From a mapping of the • 
sound field, all the essential characteristics of the source and 
in particular its strength can be estimated. If the source is 
not a point source, but has an extention in space, the far 
field distribution too can yield the directivity and strength 
distribution (this is equivalent to inverting a Green's formu- 
lation for the radiation field) . 

When the source radiates in an enclosure of finite dimensions, 
however, the situation becomes very different for two reasons. 

i) The reflected waves from the walls of the enclosure when in- 
teracting with the direct waves from the source, will interfere 
and the resultant sound field is so complicated that it becomes 
very difficult to infer from the interference pattern the 
characteristic of the source. If the signal radiated from the 
source is random, as is often the case in noise from turbulence, 
the field in the enclosure, the so-called reverberant field, is 
diffuse. If the acoustic absorption of the walls is known, then 
classic acoustic theory through the Knudsen's formula yields the 
strength of the source from the average reverberant field. But 

in this case no details can be inferred from the reverberant field. 

ii) In many instances, the back reaction of the reflected waves 

on the radiating source can cause the power radiated from the source 
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in the enclosure to be different from that in free space. In 
communication theory this is equivalent to stating that because 
the source is not a constant current source, then the power 
radiated is a function of the effective impedance of the enclosure. 

In our case, it will be assumed that the source acts as a constant 
current source, since it is reasonable to state that the sound 
field will not alter the strength of the turbulence in the jet. In 
order for the sound to interact with the turbulence (i.e., a 
third order effect in magnitude) one would need extremely intense 
sounds of an energy level comparable to the energy level of the 
turbulence. Now in subsonic flow it was seen that the efficiency 
of mechanical to acoustic energy conversion went like M 5 . Typically 
the efficiency was of the order of 10~ 3 % and under extreme condi- 
tions could approach one per cent. For supersonic flows it is 
still of the order of a per cent. Hence the interaction is ruled 
out. 

One has to point out, however, that for coherent noises (edge tones, 
etc.) the noise generation mechanism is very sensitive to feed back, 
consequently the reflected waves can be very important (because of 
phase effects) . It is not easy to generalize the case of the nar- 
row band noise produced aerodynamically and each case has to be 
studied individually. The only remark one can make, however, is 
that if the wave length is much smaller than the size of the enclo- 
sure and the walls are efficient absorbers, hopefully, we can 
approach the free field conditions. 
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In the next section on diagnostics we will describe three methods 
which will be able to allow one to obtain information on the 
acoustic properties of the source from reverberation field 
measurements . 

III. METHODS OF MEASUREMENTS AND INSTRUMENTATION TECHNIQUES 
3 . 1 IMPORTANT PARAMETERS NEEDED 

In the light of the previous discussion on aerodynamic noise the 
field parameters of interest are of two categories. The first 
refers to the source , the second to the field radiated at a dis- 
tance. The state of the theory is such that we can deduce the 
second (i.e. the noise field) from a knowledge of the first (for 
a few cases e.g. subsonic flow); but not the inverse. 

The source is described by means of: 

i) Strength of velocity fluctuations (turbulence level) 

ii) Spectrum of velocity fluctuations 

iii) Characteristic dimension for coherence (correlation length) 
and under certain conditions. 

iv) Strength of fluctuations of the fluid density 
The radiated sound field is described by 

i) Directivity (in free space) of source 

ii) Frequency spectrum for sound (pressure, or velocity in free 
space) 

iii) Acoustic power radiated. This can be inferred from (ii) 
above if measurements are in free space, if not one has to make 

use of methods described in VI to deduce the free space equivalents. 
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3.2 SURVEY OF KNOWN DIAGNOSTIC METHODS 

The measurement of the fluctuating velocity field in a turbulent 

(5) 

field has been obtained directly by hot wire anemometers and 

by velocity microphones. Unfortunately, either of these two in- 
struments have to be immersed in the flow with a consequent dis- 
turbance of the mean value of the flow. 

Several methods have been developed in which the magnitude of the 
fluctuating velocities and/or densities are measured by detectors 
which are outside the flow. Now two general procedures have been 
used according as to whether a tracer material is or is not intro- 
duced in the flow field. To prevent the tracers from separating 
from the bulk of the fluid it has to be present in extremely small 
quantities. 

3.2.1 NO TRACER-METHODS 

In these methods, the interaction of the light (e.m. wave) with the 
sound field is brought about via the modulation of the index of 
refraction (for e.m. wave) by the fluid fluctuations. This inter- 
action has been detected in a number of ways. 

a) Funk and Cikanek v ' correlate the output from two optical 
signals, originating from two lasers and which have been transmitted 
through the gas. Although their paper does not state it, the cor- 
relation they observe is directly related to the density fluctuation. 

(17) 

b) Gruber and Marhic have conducted the above experiment 

much more carefully. They also work out the theory. They show 
that the spectrum of the signal is proportional to 

(c) 2 h / e ' iW T u«> • 


■—CO 
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V q is the frequency of the e.m. wave, k = u>/U 
number corresponding to the acoustic frequency w , 
mean flow velocity and <n(x,t) n’fa' , t+x)> 

the cross correlation for the density. 


is the wave 
is the 
is 


3.2.2 TRACER-METHODS 

If a tracer gas is introduced in the field, the sensitivity of the 
diagnostic techniques can be improved. This is because the index 
of refraction (real and/or imaginary part) for the tracer is usually 
more sensitive to the fluctuations then that (i.e. index) for the 
host gas. In most cases the hot gas is air. A typical application 
of this method is described below. 

Wilson et. al., have used a crossed beam optical method to detect 
the influence of the fluctuations on the absorption of the traces 
of CO 2 that are already present in the air. 

The transmission loss for each beam is directly related to the 
local fluctuation of the density in the mechanism. Hence the cor- 
relation of the outputs from the two beams yields directly the 
cross correlation for the density. This has been denoted by T ±i in 
3.2.1(b). Wilson et. al., assume that the C0 2 is so diffused in the 
host gas, that it is reasonable to consider the density of the 
tracer to be directly proportional to that of the host gas. 

It is interesting at this point, to observe that a great deal of 
work on holographic and optical investigations in other fields seem 
to be unfamiliar in acoustics as they are rarely referred to. For 

(19) x- (20) 

instance, it was shown by Buchshaum and Gsanastein , Bekefi 
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and others that if a mild glow discharge is struck in a gas flow 

the very weakly ionized gas yields an electron gas which can be 

treated as a tracer gas. Since the ionized gas is completely 

collision dominated , there is no question but that the local fluid 

ve locity as well as the local density of the electron gas is that 

of the host gas. A great advantage of this method is that the 

chemical composition of this tracer is independent of that of the 

host gas i.e. it will always be an electron gas, unlike the case 

of Wilson et.al., where they had to concentrate on CC^, etc. Now, 

for most of these weakly ionized gases one can easily obtain 

densities of 10 11 to 10 13 /cm 3 which brings them within the reach of 

many readily available detection techniques , in particular 

( 21 ) 

i) Laser interferometric techniques: Ashby-Tephcott and 

„ ( 22 ) 

Hooper 

( 23 ) 

ii) Holographic interferometric technique: Jahodo v 

(24) 

iii) Thomson scattering ■ 

Another advantage of these methods, is that whereas in the methods 
described in 3.2.1(a) and (b) only, the correlation of the density 
is measured, the techniques of (i) to (ii) yields the density, 
while (iii) yields the velocity. 

3.2.3 LASER INTERFEROMETRIC TECHNIQUES 

The laser cavity is coupled, in series, to another cavity which 
contains the turbulent region of interest. This region is now 
outlined by a glow discharge so that the region is filled with 
diffuse electrons. A mirror at the end of turbulent region is now 
made to rock to and fro. Consequently, the cavity containing the 
turbulent region goes through resonances and anti-resonances and 
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the corresponding output from the laser is directly proportional to 
the "load". 

It is shown that, the electron density and hence that for the gas 
is proportional to the slope near the anti-resonance. By using a 
separate path (as a reference path, one can obtain the auto- 
correlation for the density. Limiting resolution times of 0.1 
|isec have been obtained with this method. 

3.2.4 HOLOGRAPHIC INTERFEROMETER 

The advantage of this method is that one can obtain a map for 
the autocorrelation of the density. Unlike 3.2.3 where the auto- 
correlation function is measured point by point. The arrangement 
for this method is shown in Fig. 1 which is self explanatory. 

(After Jahodo, Ref. 23.) 

It is known from holographic theory that in fact the intensity 
of the signal received is directly related to the auto-correlation 
density for the electrons and hence for that of the gas. 

3.2.5 INTENSITY OF OPTICAL LINE 

This method does not rely on lasers. It follows directly from the 
fact that whereas the density of electrons is proportional to that 
of host gas, so also is the density of the excited atoms (because 
the gas is collision dominated) . Consequently the intensity of the 
light emitted at a given optical wavelength is proportional to the 
density of the gas. 

When one monitors the light emitted, one measures directly a quantity 
proportional to the density. A correlation of the light output hence 
yields the correlation of the density directly as was shown by 


Bekefi et. al. 




Fig. 1. 


/ p *i \ 

Holographic Interferometric Arrangement of F. Jahodo 1 
to Map the Autocorrelation of the Density. 
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3.2.6 THOMSON SCATTERING 

Sometime ago Salpeter^^ showed that the Thomson scattering 
from the ions in an ionized gas yields a direct measurement of 
their local temperature. Since the ions, here, are cold, it is 
reasonable to assume that their fluctuating energy (mainly tur- 
bulent energy) is greater than their thermal energy. Consequent- 
ly, the scattering energy in a Thomson scattering experiment 
should be directly proportional to the auto-correlation of the ions 
velocity (and hence of that of the neutral gas) . 

This experiment, which has been very successful in plasma investi- 
gations is difficult to perform and not conceivable to be completed 
in one year. Because of the low intensity of the scattered signal 
(since the ion density is so low) one would have to use substantial 
optical energy as obtained from a Ruby laser. The experiment is 
feasible in principle and should be considered in long range 
planning. 

IV. HOLOGRAPHY AND VISUALIZATION OF SOUND FIELDS 
We discussed in section (3.2.4) a tried technique utilizing holo- 
graphic interferometry to record a map of the autocorrelation of 
the density of a gas. In this section we shall explain other 
methods of holography which transduce the sound field to an optical 
one and hence revealing visible images of the sound sources. 

4.1 BRIEF INTRODUCTION TO HOLOGRAPHY 

The principle of holography in brief is that when a wavefront 
(emitted from or scattered by an object is attended by a coherent 
background (reference wave), the resultant record of the two waves, 
known as the hologram, contains all the information about the object. 
Hence, when the hologram is re-illuminated the original object is 
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„ . _ _ . . , (29) f (30) r (31) The 

reconstructed in its true original form • me 

reconstructed image can be a scaled up or down 3-dimensional 
replica of the object, depending on the linear scaling of the 
hologram and the wavelength of the reconstructing wave relative 
to the object wave. A hologram can be made with acoustic radia- 
tion and reconstructed with an optical one and hence transducing 
the acoustic object or source to a visible one , • 


4 . 2 REFERENCE WAVES IN ACOUSTICAL HOLOGRAPHY 
The reference wave is an essential component in recording a 
hologram. Such a reference wave in acoustical holography, how- 
ever, need not be a physical one. It can be simulated electronic- 
ally ) or temporally ^ 34 h The electronically simulated reference 
hologram is made by mixing (before recording), electronically, the 
object signal with another one from a local oscillator tuned to 
the same frequency as that of the object, and properly phased to 
simulate the required wave form. The temporal reference, on the 
other hand, has a frequency that differs from that of the object 
by an amount Af. Both types are made by either a scanning de- 
tector or a properly arrayed detectors (linear or mosaic con- 
figuration) properly phased to introduce the off-set angle 0 that 
separates the object wave from the simulated reference wave, 
according to the equation sin 0 = (Af/f) • (c/v) where v a nd 

c are the scanner and sound propagation velocities. 


4 . 3 DIRECT SOUND-LIGHT INTERACTION 

A pseudo-acoustical-holography approach is based on the direct 

(35) 

interaction between coherent light and the sound field 


Such 
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interaction results in Bragg diffraction of the light wave which 
provides a real-time visualization of the sound field. The Bragg 
angle is governed by the ratio of the light to sound wavelengths . 
The reconstructed image is demagnified by the ratio of the two 
wavelengths and furthermore, its resolution is influenced by the 
size of the light source. It has been shown analytically that 
this technique is useful for visualizing single frequency sound 
sources of medium frequencies range, (outside that of aerodynamic 
noise) and hence will be considered no further . 

4.4 HOLOGRAPHIC INTERFEROMETRY 

It is well known and has been demonstrated that a single hologram 
can store several wavefronts, e.g. at different times. Hence, if 
a conventional doubly exposed, pulsed, optical hologram 
is made of the flow, its reconstruction may reveal the change 

that took place in the flow during the interval between the two 
exposures. Careful analysis of such technique in order to locate, 
and define the sound sources in the flow may prove to be quite 
beneficial . 

4.5 DIFFICULTIES OF HOLOGRAPHY AND HOW TO OVERCOME THEM 
Holograms have been made so far in a quiet environment with 
monochromatic radiation. Both of these conditions may not be 
satisfied in the 40 x 80 foot wind tunnel: the aerodynamic noise 
sources are broad band and the wind tunnel is reverberant. Further 
more, the source of noise may not be stationary. To overcome these 
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difficulties, we must use very narrow band detectors ,. well- 
focussed on a specific plane, and more or less near-instantaneous 
type recording. With such detectors one can make a temporal 
reference type hologram (the unknown signal is mixed with a known 
fixed frequency one before recording) the construction of which will 
reveal the sources of that particular band. On the other hand, if 
the detectors have a wide band characteristics, then the construc- 
tion of a temporal reference hologram is expected to reveal one of 
the following: 

a) If the noise sources are stationary and have the same single 
frequency, they will be reconstructed optically with their relative 
positions preserved. 

b) If the noise sources are stationary, monochromatic but of 
different frequencies, they will be reconstructed, with distorted 
relative positions, (equation of section 4.2). The exact relative 
positions can be determined by repeating the experiment with 
different reference waves of different frequencies, and comparing 
the various reconstructions from the various holograms. 

c) If there is only one stationary noise source having a broad 
continuous band, its image will be revealed as a defocussed image. 
The intensity distribution in the image is a function of the power 
spectrum of the source . 

d) Discrete sources of the above (c) type can be similarly 
analyzed. The degree of overlap of the reconstructed images depends 
on the band width of each source and their relative distribution in 


space . 
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e) Moving sources of narrow spectral band width produce streaks 
of constant intensity, while broad band sources produce streaks of 
variable intensity in directions other than the direction of migra- 
tion of the source. This is the most difficult of all cases and 
require further theoretical analysis to determine both the degree 
of accuracy in interpreting the results, and the experimental 
implementation . 

f) Other configurations are possible and require the analysis 
of different holograms taken with different reference waves and 
also, taken in different locations. 

The detectors in this approach can be either 

(i) A single detector 

(ii) A linear array of detectors scanning in a direction 
perpendicular to their arrangement 

(iii) A mosaic of detectors in a two dimensional lattice, or 

(iv) A series of piezoelectric crystals probed with an electron 
beam, after a fashion similar to Sokolov tube (36) . 

The choice of the appropriate configuration depends mainly on the 
stability of the sound sources. 

Another problem that may be encountered is the effect of turbu- 
lence. This, however, can be overcome by making a conventional 
type hologram with the reference wave penetrating the same tur- 
bulent flow. Experimental verification of this approach has been 

. ( 37) 

demonstrated by Aoki 
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V. ANTENNA THEORY AND ACOUSTIC NOISE MEASUREMENT 

Measurement of the coherence of the pressure variation between pairs 
of points at various vectorial separation is the tool that is use- 
ful in the study of noise radiation from extended sources in other 

( 28 ) 

subjects, for example radio astronomy . While undoubtedly the 

same statistical theory holds true in acoustics, there are differ- 
ences in practice arising from the necessary proximity of the sources 
to the measuring instruments and from the special character of the 
instruments themselves, that require attention before an experi- 
mental investigation of this tool can be recommended. 

The coherence measurements concept is independent of the hardware 
implementation, applying equally to any acoustic sensor whose output 
follows the time-varying pressure. 

From a study of this approach, we should arrive at theoretical limi- 
tations on the resolution with which noise sources can be localized 
by measurements made well outside the main flow, where the problems 
of perturbations are at a minimum. 

Data processing of information obtained by multiple probes is 
fortunately well advanced and covers a number of configurations such 
as variable-spacing and arrays of various kinds, all of which should 
be examined for possible application to, or modification for, the 
noise problem. 

VI. DISCRIMINATION AGAINST REVERBERATION NOISE 

It was mentioned in our general remarks that three methods can be 
used to discriminate against the reverberant signals. These are 
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1) The use of a pulsed flow 

2) The use of correlation techniques 

3) The use of highly directional microphones 
We discuss these briefly. 

6 . 1 THE USE OF PULSED FLOW 

If the flow is pulsed, with characteristic times less than 30 
millisecs (corresponding to time of travel of sound over the 
smallest dimension of tunnel i.e.^40 ft) , that a discriminating 
circuit gated over the same period will only receive the direct 
sound. For large jet engines this may not be practical. 

6.2 THE USE OF CORRELATION TECHNIQUES 

Correlation techniques has been used with great success in a number 
of applications. The basic idea is based on the use of two micro- 
phones. The output of these two microphones is correlated. The 
signal received from one consists of the direct signal plus signals 
which have suffered many reflections from the walls. If the signal 

to each microphone is written as 

S d (r,t) and S^r.t) 

i.e. direct and indirect, then the correlation of the two outputs is 


F(r,r 



(r,t) + S i j(r,t)][S d> 2 (r' ,t+T) 


+ s i ^(r' ’ t+ x)] dt 


Since the reverberant signals usually suffer several random reflec- 
tions and changes of phase specially if the walls act as diffuses, 
then beyond a correlation length of the order of a wave length F(r,r‘ 
has contribution only from the two coherent signals which are the 
direct signals, i.e. 
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F(.r,r' ,t) % </S c j(r,t) 


The disadvantage of this method is that it is very laborious, as the 
field is plotted point by point, and furthermore the correlation is 
made in time, which requires special long delay instruments such 
as magnetic recorder. 

6 . 3 THE USE OF HIGHLY DIRECTIONAL MICROPHONE 

The most direct method is to - discriminate against the reverberant 
noise by means of direction signals. There are two types which are 
specially recommended 

1) The tubular microphones 

2) Shaded line microphones 

The art of designing directional microphones is very old and dates 
from the decade of the 40’ s. The characteristics of such instru- 
ments are very well known. 

Tubular microphones orginally due to Mason and Marshall and 

later refined by Olson ^ ^ have a directivity which is very sensi 

tive to frequency. In fact it is expressed by 


Directivity 


1 si n 

n si n (Jin 


where ^ A ^ C0S ® ~ ^ 

" n " is the number of tubes in the bundle. L the length of the 
longest tube, 6 the angle of view and A is the wave length of 
sound. A fifty tube microphone yields a 60° lobe with a discrimina- 
tion of 30 db for L/A % 0.5. Of course the directivity improves 
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with increasing frequency. Since aerodynamic noise, however, has 
a wide band character it is important to cover a fairly large 
range of frequency. 

In the shaded line array of microphones , the line consists of micro- 
phones not necessarily equally spaced and phase and/or amplitude 
compensated. 

For the particular case considered here, the simplest array con- 
sidered is a cross array made of two lines containing m and n 

microphones each seperated by distance£d^,d 2 respectively. The 

(27) 

directivity for such an array is given by : 

sin{ (cos a) (mTrd^/A) } sin{ (cos 3) (mrd 2 /X) } 

Directivity = — — - — 

m sin{ (cos a) (ird^/A) } .n sin{ (cos 0) (TTd 2 /A) } 

a, 3 are the zenith and azimuth angle of view. 
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ACOUSTIC ANTENNA 
ELECTRONIC- STEERING. CONTROL 


MODEL 909B 


lm mmi jjm 

Model 909B electronic steering control is a unique instrument 
capable of steering two crossed linear arrays of acoustic 
transducers (microphones) around any chosen axis in space 
withount mechanically disturbing the transducers in any way. 

In other words, with two simple dials, any position of the 
plane of the acoustic antenna may be chosen. 

The maximum number of transducers which can be used with this 
model is nine, arranged in two groups (A & B) of five each. Each 
group is placed along a straight arm, and the two arms form, 
preferably, a right-angle cross. The plane of the transducers 
can be made to rotate clockwise or anticloskwise around a 
vertical axis, out of, or into its own plane around a horizontal 
axis, or tilt in any direction by the combined rotation around 
the two axes. Thus, linear, two dimensional, as well as helical 
scanning are possible. All this can be accomplished , simply, by 
appropriate positioning of two switches and properly controling 
two dials. The dials control the necessary phase shifting of 
the transducers in order to induce the eqivalence of mechanical 
rotation. The phase shifting is independent of the frequency 
of the detected radiation. Provision is made for high directivity 
as well as improved discrimination against unwanted background 
noise, such as reverbration and echoes, through Dolph-Tschebychef f 
shading. 
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SPECIFICATIONS 


Maximum number of transducers (microphones) per arm — 

Number of linear arrays 

Maximum number of transducers (total) 

Possible steering of one array independent of the other- 

Possible steering of both arrays simultaneously 

Phase shift at 1 KHz (per delay line) 

Steering angle of any linear array 

Steering dependence on radiation frequency 

Suggested optimum value of (d/A) 

Half angular resolution 

Maximum input per channel 

Power supply 


Style rack mount 19 "w x 8.75"h x 10"d (48.3 x 22.2 x 25.4 

Weight 15.25 lb (6. 

Output each linear array (arm) seperately, or 

the product of the two arrays. 


5 

2 

9 

yes 

yes 

4.7 radians 
0° to 36o° 
none 
0.75 
10 ° 

30 mv 

110 V ;60^?1 A 
cm) 

8 Kgm) 



EL-SUM CONSULTANTS 


B-3- 


AA-MODEL 909B 


DESIGN FEATURES OF MODEL 909B 
MEASURED OUTPUT 


If w (x,y) is the microphone 
strength distribution in the 
(x,y) plane, its response to 
the incident wave (Figure 1) 
is then given by a quantity 
proportional to 

/w(x,y) e ik(x coscl+y cos6) dx dy 

where k=(2TT/A),A being the 
wavelength of the incident wave. 



For a distribution of point microphones along the x & y-axes, 
the response is proportional to: 


_ ik(x cos a + y cos 8 ) 
Ew e n n n n 

n 


Ea in the x-direction, and 
n 


„ ik(x m cos a + y m cos $ ) = Eb m in the y-direction 
lw em m in m m 

m 


where n & m are the number of microphones in the x & y-directions 
respective ly . 


Model 909 (the predecessor of this model) measures (2a n ) + (Ib m ). 

This model 909B measures (Ea n ) alone, (Eb^) alone, or the product 
(Ea n )x(Eb m ) . 

Provision is made to filter <Ea n ) and (Eb m ) before their multiplication 
operation. This filtering is necessary when the incident wave 
is noisy or hetrochromatic . If the incident wave is monochromatic 
(single frequency) no filtering is necessary. 



EL-SUM CONSULTANTS 


B-4- 


AA-MODEL 909B 


STEERING OF' ANTENNA 


When the antenna is tilted by a (or (3) ,as shown in Figure (1) , 

0 •rr'FrJ 

there is a phase difference of - cos a (for (3=0, say) between 

c 

neighbouring microphones, where d=distance between microphones, 
f is the frequency and c is the velocity of sound. 

In order to steer the antenna, therefore, we shift the relative 
phase of the output of microphones by means of electronic means 
through an angle equal to — — cos a for horizontal direction 
and through an equivalent appropriate angle for the vertical 
direction. 


The block diagram for the arrangement is shown below, Figure (2) . 

The phase shifters are all pass networks designed to provide 4.7 
radians (+ or -) shifts at 1 KHz. Since the characteristics of an 
all pass network is to have constant amplitude and linear phase 
shift, the electric shift is proportional to the frequency. 

4.7 

Actually, the electric phase shift 4>= f.a , where a is an 

10 3 

adjustable number that controls the RC value and hence the phase 
shift. Now we see that the phase shift 


wd 

c 


COS Ot 


4.7 
10 3 


f.a 


then notice that the actual steering angle a (or g ) is given by 


a 


COS {- 


4.7 


d 10 



It is more convenient to use the steering direction which is 
(90 — a) or (90 - g ) . The quantity actually varies between 0 and 


and is controlled by the dials A & B . 



same 

















OBSERVER FACES ARRAY 





FIG. 3 


Phase Shift Selector Switch 


EL-SUM CONSULTANTS B-6 AA-Model 909B 



EL-SUM CONSULTANTS 


B _7_ 


AA-MODEL 909B 


The crucial point is that the above relation for a ,(or 3) is 
independent of the frequency and is completely set by the design 
parameters (absolute phase shift and d) . 

The calibration of the instrument (dial reading vs. steering 
angle) is shown in Figure (4) 


SHADING (TCHEBYCHEFF CORRECTION) 

For five transducers per arm, the weight of the coefficients 
are as indicated below. The value of d/A is 0.75 


b _ 2 

0.46 

b -l 

0.82 

b +l 

0.82 

b . 

0.46 

+ 2 

a _2 

0.46 

a -l 

0.82 

a o 

1.00 

a , 

0.82 

+1 

a 

0.46 

+ 2 
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Figure (4) 


Instrument Calibration Curve 
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CALIBRATION AND TESTING PROCEDURE 


Calibration is preferably done in an anechoic chamber. 

( 1 ) _ F 0r (d/A)=0.75 and f-500 Hz the microphones are seperated 

by d=19" (48.6 cm) , (where d is the distance between adjacent 
microphones) . 

(2) - Mount the nine microphones on a right angle cross (Figure (1)). 

Each arm of the cross consists of five microphones. Label 

the arms, "A” for the horizontal and "B" for the vertical. 

Label the microphones, a and b where n & m are 0 to 2 . 

n m 

Set the antenna to face a loud speaker 20 ft (6 m) away 
(or more) . 

(3) - Tune the loud speaker to 500 Hz and the output of any 

microphone to not more than 30 mv . 

(4) - Short circuit the filters connections at the A & B output 

terminals (found at the back) , as long as only one frequency 
is being used. 

(5) - Switch on the power supply to model 909B . 

(6) - Connect all nine microphones to their corresponding terminals 

at the back of model 909B . 

(7) - Measure the input to each microphone using its marked tap 

at the back of the instument. Adjust this input to reasonable 
value (e.g. 12 mv) , not more than 30 mv , using the adjusting 
screws above the microphone terminals T 

(8) - Turn both dials of channels A & B (located at front) anti 

clockwise until they click. The antenna is now in the (0,0) 
position with its plane parallel to the incident palne wave front. 

THE INSTRUMENT IS NOW READY FOR TESTING WITHOUT SHADING 

IF SHADING IS REQUIRED THEN INSTEAD OF STEP (7) THE INPUT 
TO THE MICROPHONES SHOULD BE WEIGHED BY TSCHEBYCHEFF 
CORRECTION FACTORS FOUND IN THE TABLE ON PAGE 7. 

(9) - Connect the output of the instrument (found at the back) 

to a chart recorder. 
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(10) - With the antenna in the (0.6) position (step 8, above) 

turn both arms A & B mechanically through 360° and record 
the output. 

(11) - Repeat step (10) while dial A is turned to a reading 

corresponding to a specific steering angle (using the 
calibration curve, mounted on the instrument, and also 
reproduced in Figure (4) of this manual. The recorded 
output should be the same as in step (10) rotated through 
the steered angle. Notice any discrepancy and make 
necessary correction. 

(12) - Repeat step (11) using dial B alone. 

THE INSTRUMENT IS NOV! READY FOR USE. 

The following steps may be made if desired. 

(13) - Check the output of a & b ; each one should be equal to 

the sum of the output n of the microphones in the corresponding 
arm. If not, check the individual microphone circuits. 

(14) - Compare the output of each arm alone with the output of 

their product, using a CRT. The product has an amplitude less 
than the amplitude of each arm alone, and a frequency which 
is double the input frequency. 

A NEW MODEL 909C provides an output frequency , same as the input. 


IT IS VERY IMPORTANT TO CONNECT APPROPRIATE VARIABLE 
FILTERS TO THE OUTPUTS OF 2a Q Zb m WHEN THE INPUT 
IS NOISY OR N0T-M0N0CHR0MATIC, THESE FILTERS MUST BE 
TUNED TO THE PROPER FREQUENCY UNDER TEST, 



